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ABSTRACT 


ÀA coupled oceanic-atmospheric boundary layer model which 
provides single-staticn prediction capability is evaluated 
relative to wind stress and oceanic mixed layer depth obser- 
vations. The model is initialized and verified using data 
obtained during the 1983 Mixed Layer Dynamics Experiment 
MIDE ee Model computation of friction velocity (u*) and 
oceanic mixed layer depth are compared with observations 
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ations. Favorable results of model u* predictions are 
achieved althougk in some cases they are slightly higher 
than observed turbulence derived values of u*. Mixed layer 
predictions are very close to the observed except in frental 
regions. Proper ABL/OBL model initialization of boundary 
layer values is critical in order to achieve favorable 


results. 
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I. INTRODUCTION 

Naval oceanographers and meterologists obtain a great 
variety cf informaticn to help tnem in making atmospheric 
and oceanic forecasts. Extensive use of large shore based 
computers, remote sensing equirment and world-wide communi- 
cations networks make the forecaster rely on remotely tased 
technology rather than local indicators and have, to a great 
degree, removed the "art" from single-station forecasting. 

However, large-scale conflict or operational considera- 
tions could deny all forecasting inputs available to the 
forecaster except local indications. Reliance on single- 
Station forecasting for naval operations during times of 
crisis or war was first addressed by Oliver and Oliver 
MIOS): 


During the last few years wartime conditions have made 
it necessary ‚for isolated combat units to issue fore- 
casts in regions where no network of meteorological 
stations could be available. Frequently the data from 
Several stations or from, reconnaissance planes are 
available, but in scme regions the forecaster must rely 
only on surface and upper-air observations made at his 
own station. This is ا‎ true in the case of 
Ships at sea. Hence, it is important to develop profi- 
ciency in extracting information from limited aerolo- 
gical data. 


While the above quote pertains to atmospheric forecasts, the 
Same could be true for ocean predictions. 

In modern warfare, isolated ships or battle groups far 
from the main force, operating in an emission control 
(EMCON) posture, must be able to predict their operating 
environment without shore suppcrt. Such predictions are 
important for the effective utilization of weapon anc sensor 
Systems and to deny the effective use of the environment to 


the opposing force. Affected systems include: navigation 
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radars, communications equipment, over the horizon (0TH) 
radars, air and surface search radars, missile and gun fire 
control radars, laser and infra-red (IR) guided nissiles, 
hull mounted sonars, towed arrays, and variable depth sonars 
(VDS). 


(ne s rv 
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Fig. 1.1 Components of the Coupled Atmospheric 
and Oceanic Boundary Layer Model. 


A coupled atmospheric and oceanic model (O'Loughlin, 
1982) has been developed to predict ocean and atmospheric 
parameters from local indications. The components of this 


coupled model are depicted in Fig. 1.1. This model is being 


12 


tested and evaluated with shipboard applications in mind. 
Accurate predictions of ocean mixed layer depths (MLD), 
thermocline jumps and gradients are essential in optimizing 
anti-submarine warfare (ASW) systems and employing acoustic 
countermeasures. These parameters are dependent on the heat 
fluxes and wind stress at the ccean surface. This thesis 
will examine how the Single-station coupled model 
predicticns of the surface fluxes and mixed layer depth 


compare with observations. 
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II. MODEL BACKGROUND ANL ANALYSIS 


ROCEDURES 


A. BOUNDARY LAYER DESCRIPTION 


Thermodynamically and dynamically forced marine atmos- 
pheric (MABL) and oceanic (OBI) boundary layers must be 
considered in coupled model description. The air-sea inter- 
face has adjacent oceanic and atmospheric turbulent mixed 
layers which effectively insulate the gquasi-geostrophic 
oceanic and atmospheric flows atove the top of the MABL and 
below the bottom or the OBL. The primary source of the 


turbulence is the velocity (shear) and buoyancy (density) 


gradients created by the exchange of mass, energy and 
momentum across the air-sea interface. Such turbulence 
generates turbulent kinetic energy (TKE) and vigorous 


vertical mixing (Davidson and Garwood, 1984). 

Under undisturbed conditions, the MABL consists ofa 
moist, well-mixed layer extending from the surface to the 
capping inversion. At the inversion, jumps in temperature 
and humidity occur with a rapid increase in temperature and 
a corresponding decrease in humidity. Within the MABL, 
equivalent potential temperature and specific humidity are 
conservative quantities. 

A consequence of the large velocity fluctuations and 


mixing in the MABL ls that energetic eddies (extending from 


the surface to the inversion) entraln warm, dry air and 
bring momentum into the mixed layer. This results in the 
upward growth of the layer. Subsidence arises from large- 


scale atmospheric forcing and tends to limit the depth of 
the mixed layer. Therefore, the change in the inversion 
height with respect to time is a function of the entrainment 


rate and subsidence. If entrainment causes the MABL to 
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extend akove the lifting condensation level (LCL), clouds or 
fog will form within the layer and greatly affect the OBL. 

In the ocean, the OBL is defined as extending from the 
surface to the top of the seasonal thermocline. Rithin the 
mixed layer, nearly homogeneous profiles of temperature and 
density are evident. While density is a function of both 
salinity and temperature, over short time scales, tempera- 
ture has a far greater effect (Miller, 1976). At the botton 
of the mixed layer, large changes ir temperature and 
salinity exist with increasing depth. This is the top of 
the thermocline region and serves as a transition from the 
mixed layer to the usually dynamically stable interior 
ocean. Idealized model profiles for both the MABL and OBL 
are depicted in Fig. 2.1. 

Radiation is ancther source of energy in both mixed 
layers. While most solar energy will penetrate the MABL, 
this is not true of the OBL. Long wave radiation is 
absorbed and emitted by the first few millimeters of sea 
water. Therefore downward turkulent heat flux is as impor- 
tant as the upward flux. Surface heat flux in the OBL is a 
function of sensible and latent heat as well as the long and 
Short wave solar radiation. This radiation is highly depen- 
dent on cloud formation in the ABL so proper ABL modeiling 
is essential for realistic ocean MLD predictions. 

While buoyancy forces create turbulence in both boundary 
layers, relatively more turbulence in the OBL is produced by 
wind stress (Davidscn and Garwcod, 1984). The MLD will 
deepen if the wind generated downward flux of momentur and 
forced convective mixing (of surface water that has been 
cooled by radiative heat loss or evaporation) provide suffi- 
cient TKE to erode the thermocline and entrain cooler water 
downward. If, on the other hard, the near surface absorp- 
tion of solar radiation produces a downward buoyancy flux 


which dampens turbulence and prcvides for net warming, the 
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Fig. 2,1 Idealized Atzospheric and Oceanic 
Boundary Layer Temperature Profile. 


MLD will shallow. Fig. 2.2 shcws the importance of mechan- 
ical mixing to the CBL and its interaction with buoyancy 
fluxes. A primary objective of this thesis will be to 
determine the effects of wind stress on the MLD. 

Several cause and effect relationships are evident. 
First, clouds can be caused by changes in the ocean surface 
temperatures which, in turn, affect the radiation budget of 
the OBL which results in surface temperature changes. 
Secondly, the stability of the MABL is influenced by air- 
ocean temperature differences. Again, cloud formation plays 
an important role. Such relationships make a coupled 
predicticn approach necessary. This coupled model uses 


iocal oceanic descriptions froma bulk oceanic mixed layer 
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Plg: 2:2 Mechanical Energy Budget for the 
Ocean Mixed Layer. 


model and the atmospheric descriptions come from an atmos- 
pheric boundary layer model. Such a coupled model has 
tactically significant implications to forecast parameters 
that are important to weapons systems. Fig. 2.3 shows the 
interrelationships between the coupled model and tactical 


products which can use model outputs (on the extrene right 
hand side). 
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B. MARINE ATMOSPHERIC BOUNDARY LAYER (MABI) MODEL 


Existing integrated MABL models are zero-order, two 
layer, mixed layer models consisting of a well mixed turbu- 
lent boundary layer underneath a relatively non-turbulent 
free atmosphere. A transitior zone (inversion) of zero 
thickness (hence the term zerc-order) separates the two 
layers. Because it has no thickness, a jump occurs at the 
inversion of the conservative parameters. The model was 
described by Stage and Businger (1981), who formulated 
entrainment energetics, and mcdified by Davidson, et al. 
(1984), who described bulk aerodynamic formulae for the 
surface layer. 


The following inputs are recuired by the model: 
(1) an initial atmospheric sounding; 
(2) the mean winds at a level within the surface layer; 
(3) the surface temperature; and 
(4) the subsidence rate. 


As formulated, ten winds can be input over the forecast 
period. The sea-surface temperature (SST) remains unchanged 
and the surface current is assumed to be zero. SST and the 
Surface current change in the coupled version due to the 
interaction with the OBL model. Inversion height, mixed 
layer values of temperature, humidity and wind are predicted 
at 30-minute intervals in both the coupled and uncoupled 
versions. The prediction sters are shown in Fig. 220% 
Procedures are the same for clear and cloudy cases except 
for entrainment computation and estimation of cloud top 
cooling. Formation of clouds cr fog, cloud top cooling and 
associated entrainment are important in the physical 


processes in the MABL. 
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Fig. 2.4 Schematic of Input, Prescription and 
Computing Steps in BL Prediction. 


The MABL model is very sensitive to the prescribed large 
scale subsidence which moderates the depth of the moist 
Marine layer and hence, cloud fcrmation. Methods which can 
be used to compute the subsidence from single station obser- 


vations are: the kinematic method; the adiabatic method; 
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and the integration of the  mcisture budget equation  (Q) 
method. Gleason (1982) found that the Q-method has the nost 
merit and was used in this thesis. However, this method 
assumes that no advection will take place. 

Integrated rate equations by Tennekes and  Dreidonks 
(1981) are used to predict the time rate of change of the 
conservative quantities and their respective jumps at the 
inversion. The vertically integrated equations fcr a 


conserved property, X, are: 
h(DX/Dt) - mix") Ce scurce (2: 1) 


h{(DaX/Dt) = hfx(2dh/0t) = (w'x') + (EO - source (2.2) 


where: 
source = o )/pC, for xX = temperature 
source = 0 for X = humidity 


and Gamma (DI) is the conservative property lapse rate above 
the inversion and F, is the net radiative flux. The 
subscripts "h" and "o" refer to irversion height and surface 
values respectively. 

The entrainment velocity parameterization is based on 
the formulation by Stage and Businger (1981) and is used to 
close the system of equations and determine the time evolu- 
tion of the inversion height. The closure assumpton is that 
the dissipation rate of the TKE is a fixed fraction (1-4) of 
the production rate. "A" is the entrainment coefficient and 
is taken to be 0.2. 

Bulk aerodynamic formulas described by Davidson, et al. 
(1984), are used to determine the surface fluxes of 


momentum, sensible heat and latent heat. They are: 


VE o (friction velocity) (2.3) 


21 


Tx 


Co 9) (temperature parameter) (249 


q* ک6‎ q) (moisture parameter) (229) 


These fluxes are given by: 





utw' = u*2 (momentum) (216) 
T'w' = uFIA (sensible heat) (2a) 
q'w' = u*g* (latent heat) (2.09 


where C, and C, are the stability dependent wind and temper- 


ature drag "m . 9 is the potential temperature and 
q is the specific humidity. The subscript "o" denotes the 
surface value. 

Because radiation is a primary factor in the OBL devel- 
opment, considerable effort has gone into calculating the 
radiation budget in the MABL. Uncertainty in background 
aerosols, distribution and cencentration of atmospheric 
absorbing gases, and cloud droplet size spectra are sources 
of error in the MABL radiation calculations. The short and 
long wave fluxes are computed separately. Short wave radia- 
tive flux is calculated using the delta-Eddington method 
(Joseph, et al.,1976) which accounts for heating of the 
mixed layer by solar radiation. Incident flux at the top of 
the mixed layer is obtained from the flux at the top of the 
atmosphere and the average transmittance in each of 15 5 
covering the spectrum from 0.2 um to 1.7un. Short wave 
extinction is due to both scattering and absorption. Water 
vapor in the atmosphere is the primary absorbing constit- 
uent. Scattered radiation formas a second short wave radia- 
tive component (diffuse solar radiation) which, in the mixed 


layer, is due to atmospheric [articles (cloud droplets and 
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aerosols). Total shcrt wave radiative flux at the surface 
is determined by computing both the direct and diffuse radi- 
ative components. Fairall, et al. (1981) gives a comprehen- 
Sive review of this method. Solar zenith angle is necessary 
in order to compute short wave radiative flux andit is 
determined from the latitude, julian day and time of day of 
the initial input. The fraction of short wave radiation 
from the sea surface, Ag, which is an important short wave 
ratiative parameter, is prescriled in the MABL model as 0.1. 

Long wave radiative flux is calculated differently for 
the cloudy and cloud free cases. In the cloudy case, the 
only clouds permitted are non-black stratus. Radiative flux 


is a function of cloud enissivity, € which depends on the 


C , 
total liquid water ccntent, W, cf the cloud. W profiles are 
nearly linear with height above the LCL (Davidson, et al., 


1984) and, along with €, can be described by: 
ol 2)89, (2.3) 
er ہے‎ exp (-a¥) (2210) 


Where p, is the density of air (1.25 X 1073 gnm/cm?), his 
the height of the mixed layer (cloud top), Z, is the LCL 
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(Mou bottom), a = 0.158 (Slingo, et al., 1982), and di is 
the liquid water content at the cloud top. 
Using the Stefan-Eoltzmann law, the net long wave cloud 


on radiation flux, I can be calculated from the cloud 


nh? 
۴۰۰۰٣6۵۶۶۶۵۲1۶27 IT, ~ and the effective radiative sky 
temperature, PS , cloud bottom net radiative flux, L,,, is 


calculated in the same way using the cloud bottom tempera- 


ture, T, , and the SSI, T,. These fluxes are given by: 


MICE. e T.) 2211) 
05 وج‎ wy 2) 
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where o is Stefan's constant (4.61 X 10711) and EE 
obtained from eguation 2.10 . The net long wave flux at the 


surface, F becomes: 


long? 


o [ > 


Pong x e, Tf Em - €) ےت‎ (<ol>) 
where T is the average temperature of the cloud. 
For the cloud free case, the net fluxes are calculated 

at the top of the mixed layer (2 = h) and at the surface (Z 
= 0) by integrating the flux emissivity profile (Fleagle and 
Businger, 1980), which is a function of the water vapor and 
temperature profiles. The net long wave flux at the surface 
for the clear sky case is: 

Pons SEC (2. 14) 
where F, and F, are the upward and downward radiative fluxes 


respectively. 


C. OCEANIC BOUNDARY LAYER (OBL) MODEL 


The OBL model is a one-dimensional, second order turtu- 
lence closure, vertically integrated (bulk) model of the 
upper ocean surface turbulent bcundary layer or mixed layer 
developed by Garwood (1977). It uses the continuity equa- 
tion for an incompressible fluid, the first law of thermody- 
nanics (heat equation), the  ccnservation of salt equation, 
an analytical equation of state, the Navier-Stokes equation 
of motion with the geostrophic component eliminated, and a 
two-component vertically integrated TKE budget. 

An understanding of the dynamics of the entrainment 
process is necessary to predict the MLD change with time. 
The stable water mass underneath the mixed layer is destabi- 
lized and eroded by TKE in the nixed layer. This TKE budget 
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is the basis for an entrainment hypothesis. A closed system: 
of equations is obtained by using the bulk buoyancy and 
momentum equations with the mean turbulent field modeling of 
the vertically integrated equations for the individual TKE 
components. 

To better define the mixing process, separate vertical 
and horizontal equations for TRE are used. Buoyancy flux 
and shear production provide energy for vertical mixing, 
With buoyancy flux being a somewnat more efficient source of 
energy for mixing due to its direct contribution to the 
vertical component of TKE. The buoyancy equation is derived 
from the heat and salt equations using a linearized equation 


of state: 
P -85[1-2«(8 - 9,) - 8(S - 5,)] (2. 15) 


Buoyancy is given by: 


nr! 


i= 9(p.- p) Zp, (2.16) 
where: 

9 = temperature 

S = salinity 

g = gravity 

p = density 


R 
li 


thermal expansicn coefficient 


B 


density coefficient for salt 


The tilde represents the total instantaneous value and the 
Subscript "o" denotes an arbitrary, but representative, 


constant value. In the short term, salinity has little 
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effect on the upper ocean density profile except at higher 
latitiudes (Miller, 1976). This leaves temperature as the 
dominant density factor. However, by using buoyancy (b) 
instead of only temperature (8) allows the model tc be 
applied in cases where evaporation and/or precipitation 
contribute significantly to the surface buoyancy flux. 


Model initialization requires: 
(1) mixed layer temperature and salinity profiles; 
(2) wind-driven horizontal current profiles. 


The initial salinity profile is made isohaline if unknown 
and the initial currents are assumed to be zero if unknown. 
lack of initial information about current and salinity 
profiles is not a serious deficiency because the model will 
evolve reasonabie transient profiles of salinity ana 
momentum after only a single diurnal cycle, and the final 
results are not sensitive tc the initial current and 
Salinity profiles (Davidson and Garwood, 1984). It is also 
possible to prescrike an upper ocean internal vertical 
velocity (upwelling and downwelling) if it is known to be 
Significant. The mixed layer depth, h, is defined as the 
shallcwest depth at which the ofserved density value, c, is 
0.02 o, units greater than the observed surface density 
value. 

At each one hour model time step, the following boundary 
conditions are required: 


e net upward turbulent heat flux at the water surface 


(sensible and latent heat fiux) plus back radiation; 


e incident solar radiation; 


e the fraction of short wave radiation absorbed in the top 


one meter of the ccean; 


e surface wind speed and direction; 
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.e clcud cover; 


e temperature (SST, dry bulk and dew-point air tempera- 


ture); and 
٭‎ precipitation (P). 


Vertical velocity at the bottor of the mixed layer nay be 
prescribed hourly. Other physical and model constants that 


need to ke prescribed are: 


"the extinction coefficient for short wave radiation 


el sorpt ion: 
e the Coriolis parameter; 


e the critical Richardson number for stability adjustment 


below the mixed layer; 
e expansion coefficient for temperature; and 
e density coefficient for salt. 


The model forecast is not particularly sensitive to these 
constants, but reasonable values have been determined by 
Gallacher, et al. (1983), and they need not be readjusted 
for geographic and seasonal variability. 


By use of bulk aerodynamic formulas, the turbulent 


fluxes of latent heat, dor and sensible heat, Oh « can be 
estimated as follows: 

E € (98E, ~- E) Uo (2517) 

Tm C (T. = T )U, (2. 18) 


Net back radiation ls estimated from the empirical equation 
(Husby and Seckel, 1978): 
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Q = 1.14 x 1077 (273.16 +1 MOE O (2. 19) 
(1 - .6C2) 


where: 


E,= saturated vapor pressure of the marine air (0.98 


corrects for salt defects) 
E,= vapor pressure of air based on dew point temperature 


T_= air temperature 


T.= sea surface temp rarür 
C = fractional cloud cover 
C= drag coefficient 


Ihe upward heat flux, Q,, is then given as: 

Qu s Ve tee ON (22208 
and the solar radiation, Qos D y 

Q, = (1- a«”)(1 - .66C3)0, (22i 


The constants "a" and "b" are adapted fron Tabata (1964) and 


the cubic cloud cover correction from Laevastu (1960). 


Clear sky radiation, Q or is given by Seckel and 7ة‎ 
(TOT) | 
Q, = A, t 4 COS 0۹+ ھ52 رھ‎ 0 ۰۶٦ [25:729 


BESIN 20 


where the coefficients (A7 AT ers) were calculated by 


harmonic representation of the values predicted in the 
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Smithsonian Meteorological Tables (List 1958) with: 
0 = (27/365) (t - 21) (2525) 


where "t" is tne julian day of the year. 

The ocean is not very transparent to solar radiation. 
In the oren ocean, approximately 50% is absorbed within the 
first meter (this fraction is called "Rf"). Absorption 
varies from region to region and is highly dependent on the 
concentration of absorbing particles such as phytoplankton, 
"Gelbstori" (yellow substance), and suspended particulate 
matter. Coastal regions absorb more radiation because of 
increased amounts of suspended farticulates. Very little of 
this absorbed radiation penetrates below the mixed layer 
because ci the low thermal conductivity of the underlying 
stable thermocline. Thus most cf this energy is transferred 
upward out of the ocean and back into tne atmosphere. What 
Short wave radiation remains peretrates the mixed layer and 
is attenuated in an exponential fashion depending on water 


turbidity. Net heat flux at the surface is given as: 
ےڈ‎ T Cu * EEO. Qe (2.24) 


Surface buoyancy (heat and salt) and momentum fluxes can 
be computed by use of the foregcing equations. Mixed layer 
turbulent temperature, salinity, velocity and Luoyancy 


fluxes are given by: 


xx cO MDC. رت‎ 
(Stw') = (P - £5, (2. 26) 
(u'w') = u*2 0 ٦ 
(But) = gie (Tut) - B(S'w') ) (285) 
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where the subscript "o" refers to the surface value. The 


friction velocity in air, ٭ن‎ ۷" +76٤ 


0 TIE Ca Vio 3 2 
už = (7 /P,) 1/2 (2.398 
where r. is the surface stress (dynes/cm2). The fluxes of 


momentum, radiation, latent and sensible heat at the sea 
surface determine shallowing (retreat) or deepening of the 
mixed layer by entrainment. If there is a positive buoyancy 
flux (Q,,,X0 and E»P) the MLD will deepen. This will be the 
case at night or whenever long wave radiative cooling plus 
the upward turbulent fluxes of heat and moisture exceed 
solar radiation. Negative buoyancy flux will result in the 
shallowing of the MLD and is caused by the domination of 
daytime solar heating (net warming) at the surface, provided 
wind stress is not large. 

A schematic of the ocean mcdel computations is shown in 
TIS LS New ocean temperature, salinity and wind-driven 
current profiles are predicted at one-hour intervals and are 
used to predict the new MLD. 


D. COUPLED BOUNDARY LAYER MODEL 


Coupling of the atmospheric and oceanic models was first 
accomplished by O'Loughlin (1982) by matching momentum, 
Sensibie heat, latent heat, and radiation at the air-sea 
interface. The atmospheric part of the coupled model 
computes the radiative, heat and momentum fluxes which are 
inputs to the ocean model. Since the ocean model nas an 
hour time step, and the atmosfheric model has a half-hour 
time step, the ocean model was inserted as a subroutine into 


the atmospheric model and is called every other atmospheric 
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Fig. 2.5 Schematic of Input, Prescription and 
Computing Steps in OBL Prediction. 


time step. Before the atmospheric fluxes can be utilized by 
the ocean model, they must be converted to the appropriate 
units for seawater. 

The ocean model prediction requires the vector surface 
stress for the ocean turbulent velocity flux, u*,?*. This 


requires that the wind fielä be broken into horizontal 
components: 
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a 
II 


iox = “(sin 9) Uo (2:31) 


U 


10y -~ (cos 9) U10 (22m 


where U, is the wind speed at 10 meters and JUjo, , Uioy are 
the east-west and north-south horizontal components respec- 
tively. The wind direction, 8, is relative to true north. 


Then u* 2 is given by: 
T - pQu*? = pauras (2223) 


the density of seawater, Pyr to be about 1 gm/fcm3. Finally, 
the horizontal components of the turbulent velocity fluxes, 


2 ;. 
5, and e y are computed Ly: 


u*2 = CoD, 2 = Cy(Uig, 2 + Uy, 2) De 25) 
ولاو“ = 2 پں فا‎ 2355) 
eu, (2. 36) 


SST is passed from the MABL mcdel back to the atmcspheric 
model for use in the next time step. Pages 2.6 shows the 
interrelationships between the atmospheric and oceanic 
models in the coupled model. 

The coupling process has the potential of significantly 
enhancing the OBL prediction qualities because of the 
improvement in the boundary conditions inherent in a coupled 
OBL-MABL systen. Longer period reasonable forecasts should 
te obtained due to the feedback between the adjacent turbu- 
lent boundary layers and the associated thermodynamic and 
dynamic adjustments in each layer. Formation of Stratus or 


fog in the MABL model is particularly important. Accurately 
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prescribed incoming solar radiation is  imperitave for 


useable CBL forecast. 


and Garwood, 1984). 
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This can only be done by 


MABL model such as has been described in Section B 


ست ی ر و ی وو ی وو ت د ی و ہے سک سو ود یک u re Be‏ ولي ل سس 


a 


use of an 


(Davidson 


E. U* ESTIMATION FROM TURBUIENT KINETIC ENERGY (TRE) 
DISSIPATION 


Direct measurement of u* on a ship is rather difficult. 
One popular technique which is used extensively at the Naval 
Postgraduate School is the dissipation method. An excellent 
review of this method as well as the bulk method is 
contained in Schacher, et al. (1981). The dissipation 
method involves the use of the turbulent kinetic energy 
(TKE) equation which, assuming steady conditions and the 
divergence term to be small relative to the  cther terms, 


becones: 
u* = ([Ekz21/[0 (2/1) - (2/1) 78 (2. 37) 


The dencminator on the right hard side of eguation 2.37 has 
Leen parameterized, where the quantity P „(z/L) = (z/L) is 
substituted by Qc following Wyngaard and Cote (1971): 


be = (14 0.5 IZ/LI 35( 2 ۶9۶۰ ٦ (2. 38) 
be = (1+ 2.5 (Z/I) 25) Ye OL 2/7 la wo (2. 39) 


Determination of the dissipation rate, €, necessary for this 
technique requires that the inertial subrange follows a 


Kolmogorov spectrum, i.e., where: 
Su (K*) = < € 26 k* -5sf/3 (2.40) 


Sy (k*) in equation 2.40 is the spectral density of the hori- 
zontal windspeed, « is a coefficient with a magnitude of 
0.52 (Champagne, et als, 1977), and k* is the wavenumber. 
Eguation 2.40 has been experimertally verified. On the R/V 


ACania a hot film sensor was used to obtain a spectrum of 
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high frequency velocity fluctuations from which the TKE 
dissipation u* was calculated. It is important to note here 
that foor velocity turbulence spectra are obtained during 
periods cf fog, rain cr drizzle. This leads to the dissipa- 
tion method being invalid. Errors may also occur if the 
relative wind is less than 3 m/s allowing condensation to 
form on the hot fiim. 

Applying Taylor's hypothesis, i.e., k* - 2xf/U, equation 
2.40 may be rearranged to become: 


€ = 240 [ S,(£f) ]%° (£ZU) s (2.41) 


where f is frequency measured in Hz. The inertial subrange 
spectral densities are integrated between two wavenumbers, k, 
and k, respectively, in equation 2.40 following Khalsa and 


Bussinger (1977): 


442 =P S,(k)dx = (3/2) « e3} (k/2%8 - k/2) (2. 42) 


Applying Taylor's hypothesis, equation 2.42 now becomes: 


g,,* 7 0.230 (eU) 73 (£74 - £775) (2. 43) 
Rearranging eyuation 2.43 and ccmbining with equations 2.37, 


2.38 and 2.39, with avon Karman constant of 0.4, we obtain: 
u* - 2.810, (Z/[U0 (zZL) J} ٭‎ [£72 - £52} ¥2 (2. 44) 
A f € 


Substituting f - 50Hz, f = 5Hz, the measurement height, Z = 
20.0n, equation 2.44 reduces to its operational form as 


follows: 


ur = 3.69 uam 


$c (z/L) ] Y? (2.45) 
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Because z/L is a function of the friction velocity, UN 
according to equations 2.38 and 2.39, the solution to equa- 


tion 2.45 involves an iterative process (Large, 1979). 
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III. SYNOPTIC SITUATION 


A. MIXED LAYER DYNAMICS EXPERIMENT (MILDEX) 


The Mixed Layer Dynamics Experiment (MILDEX) was 
conducted during the period of 24 October to 10 November, 
1983 near 34°N and 126°H. MILDEX was a multi-group and 


multi-platform experiment designed to: 


e provide magnitudes of air-sea energy exchange rates for 


use as boundary conditions in mixed layer modelling. 


e evaluate the drag coefficiert from TKE dissipation meas- 
urements and its dependence on swell amplitude ani 


direction. 


e to provide a time series, from radiosonde measurements, 


of MABL structure for model verification. 


e provide an intercomparison between radiative transfer 
algorithms developed at SEBLPES Institute of 


Oceanography and measured data. 


e evaluate the utility of SODAR as a ship-Lorne 


instrument. 


e evaluate radiative transfer models of the atmosphere 


with cioudiness as a unigue parameter. 


Meteorology data were collected on the R/V Acania as well as 
two other platforms during MILDEX. R/V Acania data collec- 
tion consisted of measurements of windspeed, temperature, 
humidity, radiation and atmospheric pressure in the surface 
layer, and sea surface temperature (measured by thermistor), 
and radiosonde derived profiles of atmospheric  windsreeds 
and directions, temperatures and humidity. Radiosonde 


launches were coordinated with satellite pass times. 
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Fig. 3.1 depicts the locations on the bow of the R/V 
Acania of mean and turbulent teasurement sensors. There 
were two levels of these measurement sensors: bow mast at 5 
meters and main mast at 20 meters. In addition, there were 
hourly observations of sea state and cloudiness. 
Conductivity, temperature and depth (CTD) casts were also 
made from the R/V Acania. Two cther vessels, the R/V Wecoma 
and R/P FLIP collected meteorolcgical and oceanic data near 
the same site as R/V Acania. During the periods of interest 
in this thesis, the R/P FLIP was always within two kilo- 
meters of the R/V Acania. The meteorological data has been 
compiled in the MILDEX report (Geernaert, et al., 1983). 


Pa AA 
q 


MET. TOWERS 


MEAN AND 
TURBULENCE 
SENSORS 











٦ 
| 


| 0 
| a 


= i 7 f M1 Y) RADIOSONDE SHED | 
anas | 

| 

| 

| 


| 
| an 
x d fT SPY e سس سما لا‎ 7 A 





|] 


len a 5 تے‎ 
Wer لپ‎ ٣ 115 7 DI) LI NENNT 5 1: 


حجہ جہ 
- 





lA. — — < - = = —— 一 -一 一 一 一 一 一 一 


— — —— 


一 —— - —  - - — ee un‏ .سے سوہ بے ہے سے ہے س لیے 


METEOROLOGICAL DATA ACOUSTICAL 
ACQUISITION CENTER SOUNDER 


(Main Deck) 





MET. TOWERS 


RADIOSONDE SHED 





Pig. “3.1 Schematic of the R/V Acania 
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Time series of wind speeds, directions, SST, air temper- 
ature, and surface pressure in Fig. 3.2 show some of the 
variability experienced during EILDEX. Frontal passages are 
evident on 29 October, 31 Octcber, and 06 November. For 
most of the experiment, the SST was greater than the air 
temperature (except cn 29 October and 09 November), so the 
surface layer was unstable or near neutral. ت0‎ 7۶۰۰٠۷09-8 
hour periods on 29 October and 09 November was tbe surface 
layer stable. 

Three periods (cases) were considered for purposes of 
this study. These cases were chosen on the basis of strong 
wind forcing of the OBL and the availability of reliable 
radiosonde data for model initialization. ۶1۰9۰29-3 2315 
the movement of the R/V Acania for all three cases. The R/V 
Acania was located near 34.0N, 125.4W during all periods. 
All times are Pacific Daylight Time unless otherwise speci- 
fied. MABL turbulent kinetic energy dissipation rate 
derived u*'s were obtained fror velocity fluctuation meas- 
urements from the R/V Acania main mast (at 20 m) and a bow 
mast (at 5 m) instruments. MLD data were obtained from CTD 
casts from the R/V Acania for Case I and by the R/P FLIP for 
Cases II and III. 

Conditions associated with these cases was dominated by 
increasing clouds, swell and winds with intermittent to 
heavy rainfall. Heavy rain shcwers occurred on 29 October 
when the first of several fronts passed through the area. 
Since the presence of fog, rain or drizzle affects the 
validity of the hot film signals used in the dissipation 
technique, care was taken to exclude affected results during 
the periods. 

Ocean conditions were also variable. MILDEX experiment 
participants reported large scale advection and upwelling 
occurred during the course of the experiment. Seni faicant 
internal wave activity was also evident and is being 


analyzed by other investigators. 
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۲۲39 3.2 Time Series of MILDEX Wind Speed and Direction, 
SST, Air Temperature, and Surface Pressure. Times are PDT. 
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Eng. 3.3 Movement of R/Y Acania: Case I (solid), 
Case II (dash), Case III (dot). 


B. CASE I (0900, 31 OCTOBER - 0900, 01 NOVEMBER) 


During this period, the R/V Acania experienced weak ccld 
frontal passage at about 1700 or 31 October. Fig. 3.4 shows 
the National Meteorological Center (NMC) surface analysis 


for Case I. The synoptic situation was that a low center 
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was located in the northeast Pacific ocean, and a high dori- 
Rated the Pacific north of Hawaii. The entire period had 
nearly 100% cloud cover and intermittent rain or drizzle. 
Several rain squalls were noted after the frontal passage. 
Clouds varied from altostratus frior to the frontal passage 
to stratus for the rest of the Case I period. No fog was 
reported, so a considerable amount of hot film data was 
available. Fig. 3.5 is a copy cf GOES WEST visual imagery. 

As shown in Fig. 3.2, winds were initially from the SSW 
with a shift to the west after the frontal passage with 
speeds varying from 7 to 9 m/s ahead of the front and from 3 
to 5 m/s behind the front. There was a steady decrease in 
surface pressure from 1018 mb tc 1014 mb. Swells increased 
from 4 to 9 ft after the frcntal passage and gradually 
decreased to 5 ft towards the end of the observational 
period. Air temperature dropped by 1.5°C after the frontal 
passage. Throughout the period, the SST was greater than 
the air temperature with the greatest difference at night as 
expected. Air-sea temperature differences varied fron -.2? 
to -2.3° C. The MABL surface layer was near neutral to 
unstable during the period. 


C. CASE II (1600, 03 NOVEMBER - 1600, 04 NOVEMBER) 


During this period, a stronger cold front than in Case I 
was approaching- Fig. 3.6 is the NMC surface analysis for 
Case II. A low center was positioned in the Gulf of Alaska, 
and a high pressure center had intensified west of Baja 
Califcrnia and south of the Acania. No frontal passage was 
observed for this case. Winds were generally from the NNW 
shifting to the west in the middle of the observational 
period and then returned to the NNW at the end while wind 
speeds fluctuated between 3 and 7 m/s. The surface pressure 


was very steady changing only by 2 mb from 1021 to 1023 mb. 
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Swells were significant and from the NW and ranged from 7 to 
9 ft. Fluctuations in air temperature never exceeded eo ee 
As in Case I, the SST was greater than the air temperature 
with the air-sea temperature difference varying from -.6?to 
-1.6 °C. In this case, the MAEL surface layer was unstable 
throughout the entire period. 

Fig. 3.7 shows GOES WEST visual imagery for this case. 
Cloud cover was heavy in the late afternoon and again at 
mid-morning on the next day. During the night, the sky was 
cloud-free. Fog, cumulus and stratus clouds were evident in 
the afternoon. Stratus with some cirrus was present the 
next morning. No precipitation was recorded during the Case 


ie period. 


D. CASE III (1600, O4 NOVEMBER —- 1600, 05 NOVEMBER) 


The coid front approaching during Case II dissipated at 
the beginning of the period and the area came under the 
influence of a weak high pressure south of the R/V Acania. 
Fig. 3.8 is the NMC surface analysis for Case III. Another 
low pressure system is to the NW with an associated cold 
Fr Ont. This second system did not appear to influence the 
weather or the seas in this case. Winds initially were fron 
the NNW and gradually shifted to the SSW as the high moved 
through the area. Speeds gererally decreased from abcut 7 
to 2 m/s, and air temperature fluctuated no more than Am 
Suríace pressure varied by 2 mt from 1022 to 1024 mb while 
sea swells were from the NW at 8 to 10 ft. As in the 
previous two cases, the SST was greater than the air temper- 
ature with the air-sea temperature difference varying from 
-1.1 to -2.1°C. Again, the MABI was unstable for the entire 
period. 

This last case had the least cloud cover due to the 


influence of the high pressure system. Fig. 3.9 shows GOES 
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WEST visual imagery for Case III. At the beginning of the 
period, 50% cloud cover, consisting mainly of stratocunulus, 
was evident. This changed to cirrus and finally a clear 
hight-time sky. On the morning of 5 Nov. there was 30% to 
50% cloud cover (mostly cirrus and some cumulus) with 
clearing in the afternoon. Nc precipitation was recorded 


for this period which was good for hot film performance. 
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IV. DATA ANI RESULTS 


A. APPROACH TO FORCING U* INTERPRETATION 


For each of the three cases, the u* associated with the 
model predictions was compared with the TKE dissipation rate 
(€) derived u* andthe model MLD was compared with the 
observed MLD. As described previously, the epsilon (€) 
estimates necessary fcr u* calculations by this method were 
obtained from hot film measurements on the bow (5 m) and 
TMT 03و‎ 5.59 (20 m) of the R/V Acania (Fig. 3.1). 

The MABL model is configured for ten input winds, so the 
ten greatest and least winds were chosen from each 24-hour 
period. These will be referred to in the discussion as 
"peak" winds. Linear interpoiation was used for data pcints 
between these ten pcints. However, a forecaster at sea 
would not predict high and low winds, but would most likely 
forecast average winds over the period. Therefore, another 
model run utilizing 2.5 hour wind averaging was input over 
the same periods to see if the use of averaged winds would 
have any deleterious effect on the MLD computations. 


The purpose of these comparisons is to evaluate atmos- 


pheric forcing (u*) in relaticn to predicted and observed 
Pots. There are three u*'s used in this study: those 
estimated by bulk fcrmulations using peak winds, by bulk 


formulations using averaged winds, and by dissipation esti- 
mates. In the following discussion, the first two sets will 
be referred to as model u* (peak or averaged) and the last 


set as dissipation u*. 
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B. CASE I (0900, 31 OCTOBER - 0900, 01 NOVEMBER) 


This first case resulted in the poorest comparison of 
predicted and observed values of any of tne three cases. 
Initialization data is Summarized in Table I. Table II 
lists the 2.5 hour averaged winds that were used for subseg- 
uent runs. In Fig. 4.1, a comparison of the two wind inputs 
is shown. Fig. 4.2 shows a ccmparison of u*'s. ۸ھ‎ 1177 
culty in comparing model and instrument derived values for 
this case was instrument perfcrmance and equipment breax- 


down. Both of these problems were encountered. 


| a ee 
| TABLE I | 
| Case One Initialization | 
| 
Atmospheric Values | 
Bee anes )ٰ ۳ sea Temp (C): 18.64 
D 17: 9 Spec. Hum: 10.47 
i28. EM Lvl: J351. 74m Advection: 0. 
Ocean Values 
Mixed Layer Depth: 34. 5m Temp. Jump (C): -.16 | 
Below Layer Gradient: -.0155 
| Winds | 
| Time (PDT) speed (m/s) Direction 
0900 823 181.0 | 
| 1120 95 176.0 | 
1345 7.4 2050 
1745 Sen 1999 
2030 3.4 213.0 
2130 6.4 215.0 
0140 35 259.0 | 
| 0345 6-1 247.0 
0505 2.4 2 309 | 
| 0720 55 294.0 | 


Rain, drizzle, and low relative wind speeds affected the 


accuracy of both the bow and main mast dissipation derived 
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Case One Averaged Winds 
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Locat Time (Hours) 





ime (PDT) Speed (m/s) 8 ٤ 
0900 8.3 181.0 
1130 8.9 178.0 
1400 8.2 200 
1630 8.6 196.0 
1900 4.6 263.0 
2130 5.4 321.20 
0000 5.7 253.0 
0230 4.5 253.0 
0500 14 2 De 
0730 4.9 295.0 





en Comparison of Case I 
ni 


tial and Averaged Winds. 
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Fig. +2 Case I Priction Velocities 


u* values. Post-processing remcval of the bad data resulted 
in significant data gaps. For the bow mast, dissipation u* 
values are missing from 1616 tc 2315 and from 0604 onward. 
Only 12 useable observed data points were available (as 
compared to 48 for the model) fcr the period. The field log 
indicates a lack of confidence in main mast spectra by the 
observers onboard due to Spurious voltage readings and 
connector problems, so data from it were not included at all 
in this comparison. 

Fig. 4.2 shows the peak wind and averaged wind model u* 
and the dissipation u* values. As expected, the model u* 
values are a function of the differences in the input winds. 
Where the averaged winds are higher, the u*'s are higher. 


The dissipation u* values obtained from the bow mast have 
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the same trend but are lower than both model u*'s in the 
pre-frontal period (before 1700). They are close to the 
"peaked" wind model predictions for the rest of the time. 
All of the u* dissipation values are lower than the averaged 
wind model predictions. 

These results of model and dissipation u* values were 
not completely expected for the frontal regime. Previous 
comparisons of bulk and dissipation u*'s in frontal regions 
(Geernaert, et al., 1985) indicated that the dissipation 
u*'s were higher than bulk values prior to a frontal passage 
and then the magnitudes reversed after the passage. Here 
the dissipation u*'s remained lower tnan the bulk value 
throughout. Generally, the study is revealing that small 
pre-frontal wind speeds (less than 6 m/s) produce greater u* 
deviations (turbulent vs bulk). Additionaily, if the swell 
and wind oppose each other, greater u* deviations were also 
evident. These were minimized if the wind and swell were at 
right angles to each other. The reasons for this are not 
entirely clear. However, it has been postuiated (Geenaert, 
et ai., 1985) that the directicn and magnitude of the swell 
relative to the wind (primarily the interaction of wind 
stress with short wave elements on long wave crests) affects 
the roughness length, Zoe which in turn affects u* and Cae 

No attempt will be made with these results to examine 
how the swell interacted with the wind and how this could 
have modified the u* values. Perhaps this case did not have 
the same kind of swell-wind interaction as in previous 
studies. This is an area for further investigation. 

Fig. 4.3 compares the model and observed  MLD's. The 
model prediction was based on kulk u* values shown in Fig. 
ur. In this case, model results (solid line) early in the 
period do not agree with the observed MLD (dashed line). 
While both curves nearly coincide at the start, the model 


curve shallows much more quickly although it reaches’ the 
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proper depth. Then the model deepens rapidly and stays 


Significantly deeper for the rest of the case. An MLD error 
of up to 10 meters can be tactically significant for naval 
forces in predicting sonar perfcrmance. 

The early Shallowing may not constitute a model 
prediction error. A CTD cast would give the thermal MLD 
while the model predicts the turbulent boundary layer depth. 


It takes a finite amcunt of tine for the temperature field 


to adjust to the turbulent field. The expression for the 
124 7٦٥ 

L |a hi o OS (8.1) 
where: 

T = temperature 
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t - lag time 
Q,= surface heating 


Thus based on equation 4.1, the greater the u*, the greater 
the lag and the greater the heating (0,), the less the lag. 

Using 2.5 hour averaged winds showed no improvement in 
MLD prediction (Fig. 4.3 dotted line) and-in fact gave a 
poorer prediction. While the two wind inputs (Fig. 4.1) are 
quite different, the output MLI's are quite close together 
except at 1300. This is a pre-f£rontal period of strong 
winds and, from 1200 to 1700, the averaged wind speeds are 
higher than the initial winds. Since MLD is a function of 
the wind speed cubed, the strcnger winds would give deeper 
MIXING This is again true for the period fron 2300 to the 
end of the run where the averaged winds are generally higher 
than the initial winds. Overall, one can conclude that for 
this case both the "peaked" and averaged winds give poor MLD 
predictions and these could have adverse effects for naval 
operations. However, if the dissipation u* values for 0900 
to 1700 (which are lower than the bulk values) had been used 
by the model, the mcdel MLD wculd have been closer to the 
observed. 


C. CASE II (1600, 03 NOVEMBER - 1600, 04 NOVEMBER) 


The coupled model MLD prediction agreed well with the 
observed MLD changes in this case. The modei predicted 
unstable ABL conditions throughout the period which agreed 
well with observations. Table III contains the initializa- 
tion data for this case and the 2.5 hour averaged winds are 
listed in Table IV. Fig. 4.4 compares the two input winds 
for this case. 

This was a period with unreliable dissipation u* data. 


Eg. 4.5 is a plot of the various friction velocities. 
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TABLE 1II 


Case Two Initialization 
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Winds 
Time (PDT) Speed (m/s) Direction 

1600 4.0 040.0 
1745 6 310.0 
2120 8.6 036.0 
2230 4.7 31920 
2300 6.4 39320) 
0135 2.4 20 
0000 10 281-0 
0705 340 350.0 
1005 Tae 008.0 
1420 4.6 336. 0 
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While fog or rain were not a factor in dissipation measure- 
ments, equipment failure was a problen. This led toa 
reduction of dissipation u* data. The bow mast data were 
not available from 1700 to 2300, 2330 to 0217, and 0900 to 
1600 when a major equipment repair was accomplished. Bow 
dissipation u* data from 12 thirty minute periods were 
available. Main mast dissipation measurements were not 
available from 1600 tc 0752 and from 0822 to 1340. It was 
also absent after 1410 and only data from 4 thirty minute 
periods were available. 

At the beginning of the period, the bow mast dissipation 
u*'s are close to the bulk derived (model) values. However, 
during the middle part of the period, the bow dissipation 


values are consistently lower than the bulk (model) although 
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TABLE IV 
| Case Two Averaged Winds 
0 2) Speed (m/s) Direction 
1600 329 040.0 
1830 4.5 0د‎ 0 
2100 5. 7 IO 
2399 4.3 308.0 
0200 36 274.0 
0430 61 282.0 
0700 4.6 307:0 
0930 6.6 dog) 
1200 6-6 348.0 
14 30 5.4 341.0 | 
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Fig. 4,5 Case II Friction Velocity 


they tend to have the same time variations. The two bow u* 
dissipation values near 0800 are close to the model values. 
It is apparent from Fig. 4.6 that the model predicted 
the MLD well enough for Naval tactical applications. The 
observed MLD (dashed line) shows sone differences from 0200 
to 1000 which are probably due to internal wave activity, 
advection or aS a result of the ship moving to a different 
area. Observed MLD data was not available for 1100 to 1300. 
This contributes to the abrupt MLD rise at 1400. The model 
(solid line) predicted the MLD guite well especially for the 
deepening part of the period. Better agreement between 
observed and model shallowing would probably be evident if 


MLD data would have been availatle for 1100 to 1300. 
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Fig. 4.6 Case II Mixed Layer Depth 


Comparison of the peaked winds with the 2.5 hour aver- 
aged wind MLD's yields some interesting results. During the 
deepening phase, the averaged winds consistently give a 
Shallcwer MLD (dotted line) than the initial winds. At this 
time the "peaked" winds are consistently stronger (Fig. 4.4) 
than the averaged winds (except at 0100) and would produce 
more vigorouS mixing. However, during the shallowing 
portion, the averaged wind prediction is significantly 
deeper than the peaked wind prediction. This time the aver- 
aged winds are predominantly stronger than the peaked winds. 
Again the predicted MLD results are expected since the MLD 
is proportional to  u*?. If the dissipation u* values for 
0200 - 0600 would have been used by the model, the predicted 


MLD would have been shallower fcr that period and closer to 
the observed. 
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D. CASE III (1600, 04 NOVEMBER - 1600, 05 NOVEMBER) 


As in the previous case, the coupled model MLD 
predictions agreed well with the observed conditions. ABL 
conditions were unstable throughout the case. Table V lists 
the initialization data utilized for this model run, and 
Table VI lists the 2.5 hour averaged winds. Fig. 4.7 


compares the two input winds. 


— I 


| 


| 

| TABLE V | 

| Case Three Initialization | 

| Atnospheric Values | 

Pressure: 1022.3 mb sea Temp Cr 18270 | 

] PG Mes 17. Spec. Hum: . 48 | 

Lift. Cond; LVI: 508 0m Advection: 0. | 

Ocean Values 
Mixed Layer Depth: 11.5m Temp; "Jump ٦ | 
Below Layer Grádient: -.0053 

. | 

Winds | 

| Time (PDT) Speed (m/s) Direction | 

| 1600 52 340.0 ] 
| 12755 6.5 338.0 
2100 5.6 35120 
| 2209 گے‎ EO 
0046 du 7 338.0 

| 0155 E 32 0 | 
0445 2.6 020.0 

| 0715 2-1 296.0 | 
0915 4.8 333.0 

1045 24 271.0 | 





This case had no fog related problems affecting the 
dissipation data but for some periods, the low relative wind 
speeds degraded dissipation reliability. Several large 
dissipation data gaps resulted. For the bow mast, data gaps 
exist for 1600 to 0018, 0100 to 0375, ١+ zon 
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TABLE VI 


Case Three Averaged Winds 


me (PDT Speed (m/s Direction 
222629 2222 | 2 5( = Sot. EE 
1830 6.0 340.0 
2100 5.6 351-0 
2330 4.3 346.0 
0200 Sp 354.0 
0430 247 004.0 
0700 28 311.0 
0930 J 2 272920 
1200 2 5 264.0 
1430 DO 235-0 
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Fig. 4.7 Comparison of Case III 
Initial and Averaged Winds. 
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Fig. 4.8 Case III Friction Velocities 


1331 onward. Only seven reliable dissipation values were 
calculated for the bow mast. Data was not available for the 
Main mast from 1600 to 0018, 0100 to 0815, and from 1410 to 
the end of the period. This resulted in four thirty-minute 
dissipation averages. Fig. 4.8 compares the different fric- 
tion velocities. 

For the entire period, the peak and averaged wind model 
u*'s are very close due to the small wind differences 
between the two (Fig. 4.7). Both the main mast and bow mast 
dissipation u*'s agree reasonably with  u*'s from the peaked 
wind model run. This is particularly true in the second half 
of the period. Some variation is expected since the dissi- 
pation u* comes from instantaneous measurements, while the 


model uses interpolated values Letween the input winds. 
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Fig. 4.9 Case III Mixed Layer Depth 


The model also predicted the MLD's well (Fig. 4.9) 
although there iS greater deviation from the observed 
(dashed line) than in Case II. As before, the actual MLD 
ShowS some internal wave or advective features (particularly 
around 0000 and 0400). Good observed MLD data was available 
for the entire period. 

The initial model run (solid line) agrees well in the 
initial stages but deviates during the above mentioned 
internal wave or advective events. The sharp "dip" at about 
0930, not present in the model cutput, is an artifact intro- 
duced by use of a spline interpolation in the plotting 
routine. The model also shalloss to the proper depth but it 


occurs about one hour early. 
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Ccmparing the averaged winds MLD with the peak winds MLD 
shows that both are very similar. The plot of input winds 
(Fig. 4.7) shows that both the initial and averaged winds 
were very close throughout the period except between 0800 
and 1000, so both model runs should produce similar MID's. 
After 1100, the averaged wind MLD (Fig. 4.9 dotted line) 
Shows a double maxima during shallowing. It is unknown at 
this time why the model predicted such maxima. They have 
been slightly exaggerated by use of the spline interpolation 
but are present in the model output. As with the initial 
wind MLD, an anomolous "dip" is observed at 1030 also due to 
the spline interpolation. Both the initial and averaged 
wind MLD's are within operatioral tolerances for Navy use 
but the forecaster should be well versed in the model prop- 


erties and limitations. 
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A. CCNCLUSIONS 


The coupled model has been subject to various past anal- 
yses and performed well. This study has examined the 
performance of the model in regimes that stretch its capa- 


bilities. Several important conclusions can be drawn: 


(1) The model overpredicts clouds which affects short and 
long wave radiative effects on boundary layer 


dynanics. 


(2) Good comparisons of bulk and dissipation u* for 


unstable conditions are made. 


(3) In near neutral or slightly unstable regimes, u* 
predictions vary systematically probably due to swell 


ana sandzstress interaction; 


(4) Better predictions are achieved if peak winds instead 
of averaged winds are used. However, a forecaster 
would only have average wind prediction estimates 
available at the time of operationally initializing 
the model. 


(5) Given reasonable fluxes cf short and long wave radia- 
tion, the model predicts routinely observed varia- 


tions in ocean mixed layer depth. 


(cet 1S dê potentially pewerful tool and an aid in 
single-station forecasting to a Naval Oceanographer 


who is aware cf the model limitations. 


The model performed very well within its limitations. 
However, further work needs to be done before it can be 


implemented in an operational rcle. 
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B. RECCMMENDATIONS 


In order to facilitate operational usage, several areas 


should be explored. Some of the aforementioned prorlems 
could re alleviated by attempting the following 
recommendations: | 


(1) Continue research and development of the ا٥۵‎ 
model with direction toward operational 


implementation. 


(2) When clouds forn, the model assumes them to extent 
from horizon tc horizon. This is unrealistic and a 
statistical probability of percentage cloud cover 
should be added. Partial cloud cover would also 
alter short and long wave effects on poundary layer 


dynamics. 


(3) Incorporate swell and wind interaction effects into 


u* calculations based on bulk formulae. 


(4) Continue to analyze MILLEX data, particularly for 
periods where there are more complete dissipation 


u*'s. 


(5) Plan for model coupling to other predictive aodels 
such as IREPS and ICAPS. 


(6) Involve students in at-sea data gathering experi- 
ments. This would enhance understanding of thesis 
research and minimize "false starts" due to 


attempting cases with faulty or limited data. 


The model has tremendous potential. Performing further 
analysis and tweaking the model will result in a powerful 
tool for the fleet; cne that can be used in times of crisis 
when extensive satellite networks may not be available for 


use by forecasters. 
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